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Abstract 

We explore numerically, analytically, and experimentally the relationship between quasi-normal 
modes (QNMs) and transmission resonance (TR) peaks in the transmission spectrum of one¬ 
dimensional (ID) and quasi-lD open disordered systems. It is shown that for weak disorder there 
exist two types of the eigenstates: ordinary QNMs which are associated with a TR, and hidden 
QNMs which do not exhibit peaks in transmission or within the sample. The distinctive feature 
of the hidden modes is that unlike ordinary ones, their lifetimes remain constant in a wide range 
of the strength of disorder. In this range, the averaged ratio of the number of transmission peaks 
Wes to the number of QNMs A^modi Wes/-^mod) is insensitive to the type and degree of disorder 
and is close to the value ^2/5, which we derive analytically in the weak-scattering approximation. 
The physical nature of the hidden modes is illustrated in simple examples with a few scatterers. 
The analogy between ordinary and hidden QNMs and the segregation of superradiant states and 
trapped modes is discussed. When the coupling to the environment is tuned by an external edge 
reflectors, the superradiace transition is reproduced. Hidden modes have been also found in mi¬ 
crowave measurements in quasi-ID open disordered samples. The microwave measurements and 
modal analysis of transmission in the crossover to localization in quasi-ID systems give a ratio of 
Wes/Wnod close to In diffusive quasi-ID samples, however, Wes/Wnod falls as the effective 

number of transmission eigenchannels M increases. Once A^mod is divided by M, however, the ratio 
Wes/fVmod is close to the ratio found in ID. 
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I. INTRODUCTION 


Two powerful perspectives have helped clarify the nature of wave propagation in open 
random systems. One of them, relates to the leakage of waves through the boundaries of 
the system and can be described in terms of quasi-normal modes (QNMs), which are the 
extension to open structures of the notion of normal modes in closed systems [l-^. The 
eigenfrequencies of the QNMs are complex, with imaginary parts that are the inverses of 
the lifetimes of the QNMs. The second perspective is that of transmission through random 


systems 
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12j. For multichannel samples, transmission is most conveniently described 


m 


terms of the transmission matrix, t, whose elements are held transmission coefficients 
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15[ |. The transmittance is the sum of eigenvalues of the Hermitian matrix . Some of these 
eigenvalues are close to unity even in weakly-transmitting samples [3,Q,H,[l3|- Knowledge 
of the transmission matrix makes it possible to manipulate the incident wavefront to enhance 
or suppress total transmission through random media 18l-l22l| and to focus transmitted 


radiation at selected points 23|]. The control over transmitted radiation can be exploited to 
improve images washed out by random scattering and to facilitate the detection and location 
of objects (23|. The great potential of such algorithms for a host of practical applications 
has recently attracted attention in both the physics 231 and mathematics communities |l| 
and references therein. 

In open regular homogeneous systems (e.g. single quantum potential wells, optical cavi¬ 
ties, or microwave resonators) each peak in transmission, or transmission resonance (TR), is 
associated with a QNM ( 2^ and references therein), so that the resonant frequency is close 
to the real part of the corresponding eigenvalue. However, despite extensive research and 
much recent progress the connection between QNMs and TRs in disordered open systems 
still requires a better physical understanding and mathematical justification. 

To this end, it is instructive to look for insights in ID systems. It is well-known 3, ^ 
that the transmission of a long enough ID disordered system is typically exponentially small. 


At the same time, there exists a set of frequencies at which the transmission coe 
a local maximum (peak in transmission), and some of these are close to unity 


Scien 
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has 


27|. In 


ID, each peak is associated with an eigenstate which is a solution of the wave equation with 
outgoing boundary conditions (a pole of the S-matrix). 

Quite surprisingly, much less are studied the properties of QNMs in ID systems with 
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weak disorder where the localization length is smaller than the size of the sample. In this 
paper we show for the first time that in completely open one-dimensional disordered sys¬ 
tems, two different types of QNMs can exist: ordinary QNMs, associated with resonant 
transmission peaks and hidden QNMs unrelated to any maxima in the transmission spec¬ 
trum. The hidden modes exist due to random scattering and arise as soon as an arbitrarily 
small disorder is introduced. The imaginary parts of the eigenfrequencies of hidden QNMs 
vary with increasing disorder in an unusual manner. Typically, stronger disorder leads to 
stronger localization of modes with eigenfrequencies that approach the real axis. However, 
the imaginary part of a hidden mode's eigenfrequency, depending on the boundary condi¬ 
tions, either is independent of strength of random scattering or even increases from the onset 
of disorder. Surprisingly, the average ratio of the number of ordinary modes to the total 
number of QNMs in a given frequency interval is independent of the type of disorder and 
remains close to the constant -^2/5 over wide ranges of the strength of disorder and of the 
total length of the system. The value \/ 2/5 follows from the general statistical properties 


of random trigonometric polynomials 28[. As the scattering strength and/or the length of 


the system increase, hidden QNMs eventually become ordinary. 

The situation is different in multi-channel random systems in which a genuine diffusive 
regime exists. The degree of spectral overlap is expressed in the Thouless number, <5, which 
gives the ratio of the typical width dp and spacing Au of QNMs, 6 = du/ Au jsl, Isl, The 
typical linewidth is essentially equal to the field correlation frequency over which there is 
typically a single peak in the transmission spectrum. The density of peaks is therefore l/Su. 
On the other hand, the inverse level spacing l/Au is equal to the density of states (DOS) 
of the medium. Thus the ratio A^res/Amod can be expected to be close to Au/Su = 1/5 


for diffusive waves. The localization threshold lies at <5 = 1 
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12l |: 5 may be much 


larger than 1 for diffusive waves so that Aj-es/Amod ~ 1/(5 and may be small. For localized 
waves, the number of channels that contribute effectively to transmission, M, approached 


unity and transport becomes effectively one-dimensional 


29l |. For example, the statistics 


of transmittance are then in accord with the single parameter scaling hypothesis SOj]. It 
is worth noticing that although the statistics of the eigenstates of disordered systems is a 


subject of intensive investigations for already more that two decades (see, for example M 
^), the statistics of the transmission resonances (peaks in transmission spectra) is much 
less studied. The comparison of these two is a challenging problem for future investigations. 
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Here we find that a connection can be made between the present ID calculations of A^res/-^mod 
and measurements in multichannel diffusive systems. This is done by comparing ratio of 
A^res to the number of QNMs divided by M, MiVreg/A^mod in multichannel systems to the 
ratio iVres/fVmod in ID, where M = 1. 

II. QUASI-NORMAL MODES OF OPEN SYSTEMS 

We hrst consider a generic ID system composed of iV + 1 scatterers separated by N 
intervals and attached to two semi-inhnite leads. The eigenfunctions 'ijjm{x,t) are solutions 
of the wave equation satisfying the outgoing boundary conditions, which means that there are 
no right/left-propagating waves in the left/right lead. Each eigenfunction is a superposition 
of two counter-propagating monochromatic waves The eigenfunction in the 

jth layer, is equal to and the amplitudes in adjacent layers are 

connected by a transfer matrix. The wave numbers km are complex-valued and form the 
discrete set = k'^ — ik'^, k'^ > 0, so that the frequencies = ck^°^\ These 

eigenfunctions are QNMs. Note that all distances hereafter are measured in optical lengths. 

In what follows, the scatterers and the distances between them are characterized by the 
reflection coefficients rj = ro + Srj and thicknesses dj = do + Sdj, respectively. The random 
values 6rj and ddj are distributed in certain intervals, and ^ ddj = 0. The last condition 
means that the length L of the sample is equal to Ndo- 

To explicitly introduce a variable strength s of disorder, we replace all reflection coeffi¬ 
cients by STj, and assume (unless otherwise specihed) that the coefficients Vj are homoge¬ 
neously distributed in the interval (—1,1). This enables to keep track of the evolution of the 
QNM eigenvalues kln°^^ as the disorder increases. The condition ^ Sdj = 0 ensures that any 
random realization with the same N contains the same number of QNMs N^od = ^kL/ir 
in a given interval Ak of the wavenumbers. 

At the beginning, let us consider the QNMs of a regular resonator of the length L = Ndo 
assuming that all reflection coefficients except ri = and r^r+i = are equal to zero. In 
this case the real and imaginary parts of the QNM eigenvalues are 



when vlTr > 0, 
when tlTr < 0, 


( 1 ) 
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FIG. 1. Spatial distribution of the intensity Im{x) of quasi-normal mode (thin light gray curve) and 
Im{x) (thick red curve) in a regular resonator with (a) symmetric (|ri| = |r/j|) and (b) asymmetric 
(kil < kfil) walls. 




( 2 ) 


where m = 0,1, 2,... 

In what follows, instead of the intensity of the mth mode, Im{x) = 'ipm\x) + ^|Jin \x) 
we consider the quantity Im{x) = (x)p + (x)p, which is Im{x) averaged over fast 

oscillations caused by the interference of the left- and right-propagating waves. Examples 
of these functions for resonators are shown in Fig. [T^,b. There In{x) is distributed along 
the system as Im{x) oc cosh[2/c"(x — x*)], where x* = L[1 — ln(|r/j/rL|)/ln(|r/jrL|)]/2. When 
\xl\ = \xr\, the minimum of the intensity is located at the center of the system, and in an 
asymmetric case shifts to the boundary with higher reflection coefficient. This property will 
be used when analyzing the properties of the QNMs of the disordered system. 

In a disordered sample, the reflection coefficients r* are random and scaled by the pa¬ 
rameter s. The evolution of the eigenvalues /c(™°‘^)(s) as the parameter s grows shows that 
QNMs separate into two essentially different types. There are ordinary QNMs whose life¬ 
times, dehned by the value of 1/fc", increase monotonically with s. Simultaneously, there 
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are “hidden” QNMs (the origin of this term will be explained in the next section), whose 
lifetimes are substantially smaller than the lifetimes of ordinary QNMs and remain constant 
when s varies over many orders of magnitude. Figures [2] and [3] show trajectories of the 
QNMs’ eigenvalues = k'{s) — ik"{s) as the parameter s grows, and dependencies 

k"{s). 



FIG. 2. Motion of QNMs’ eigenvalues under disorder strength growth. Eigenvalues move from 
above. Some eigenvalues (examples are marked by “a”, “b”, and“c”) stop their motion and “stand 
still”. 

Our numerical calculations show that when external reflectors are added at the edges of 
the sample, the imaginary parts of some of the hidden modes increase with the strength of 
disorder. 

The spatial distributions of the intensity I{x) along the system are also different for 
ordinary and hidden QNMs. The evolution of I{x) as the strength of the disorder s grows 
is shown in Fig. 01 

Initially, when s is so small that values of k" are almost equal for both types of modes, 
the distributions I{x) are practically identical and evolve in the same manner: the minimum 
is placed near the center of the sample, and slopes (which are k") decrease as the disorder 
strength s grows (see Fig. 0]). These distributions are similar to the distribution of the 
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FIG. 3. Variation of k"{s). Some k” (marked by the same letters as in Fig[2|) are independent of 
s in a broad range (several order of magnitude!) of s variation. Note that length of the system 
is equal or larger than localization length when s > 0.1. Red line shows the dependence k''{s), 
described by Eq. (1^ 

intensity in the regular resonator with a small imbalance between the reflection coefficients 
and of the resonator walls. 

When k''{s) of the hidden mode “a” reaches its plateau (see Fig. |3]), the minimum of 
its distribution shifts from the center, as in the resonator with strong imbalance between 
the reflection coefficients and rn. The slope of the distribution of hidden modes remains 
constant {k” is independent of s on the plateau), whereas the slopes of all ordinary modes 
are equal and continue to decrease as the parameter s grows. The difference between the 
distributions of ordinary and hidden QNMs is that the ordinary modes are concentrated 
near both edges of the system, while the hidden mode is nestled at one edge. 

It is important to stress that this separation of the QNMs into two types occurs when 
the disorder strength s is small so that the localization length ^\oc is large relative to the 
system length L, i\Qc 3> L. Thus, this phenomenon is not related to Anderson localization, 
but, as it will be shown below, manifests itself also when £\oc <C L. 

Notwithstanding that at s —)■ 1, the lifetimes of all hidden modes increase, these modes 
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FIG. 4. Intensity distribution I(x) along the system for ordinary QNMs “ol” (a) and “oij” (b), 
and hidden QNM “a” (c) for different values of the disorder strength s. Notations for the modes 
are the same as in Figs. [2]and[3l 


are much more resistant to disorder: they become localized at far stronger disorder than 
ordinary states. As can be seen, for example, in Fig.|3l at s ~ 0.4 —0.5, when L ~ 15 —20£ioc, 
the difference between the imaginary parts of ordinary and some of hidden QNMs is about 
of one order of magnitude. 
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III. TRANSMISSION RESONANCES IN ID SYSTEMS 


We now consider the transmission of an incident wave through the system. The wavenum¬ 
bers at which the transmission coefficient reaches its local maximum and the corresponding 
helds inside the system are transmission resonances (TR). The QNMs and TRs are interre¬ 
lated. In what follows, we explore the relation between QNMs and TRs, in particular, study 
the differences between the spectra of TRs and QNMs. 

It is easy to show that in a resonator, the wave numbers of the transmission res¬ 
onances coincide with the real parts k'^ given by Eq. o. and there is a one-to-one corre¬ 
spondence between QNMs and TRs so that the number of resonances Wes is equal to the 
number of QNMs, Wnod, in a given frequency interval. The same relation also exists in 
periodic systems (periodic sets rj and dj) 


In disordered systems, the relation between QNMs and TRs is quite different. While each 
TR has its partner among the QNMs, the reverse is not true: there are hidden QNMs that 
are not associated with any maximum in transmission as shown in Figs. |5]and[6] 



FIG. 5. Transmission spectrum T{k) at s = 0.1 (L ~ Wc)- The black solid (dashed) vertical 
lines indicate the k'^ values of the hidden (ordinary) QNMs. Every maximum in the transmission 
spectrum can be associated with ordinary QNM 2,3,4), whereas hidden QNMs ( 77 ^ 1 , 5 , 6 ) are not 
associated with the transmission resonances. 

Figure [7] illustrates another fundamental difference between the ordinary and hidden 
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FIG. 6. Variation of the wave numbers (red crosses) and k'{s) (closed blue circles) with 

the strength s of the disorder. QNMs are numbered as in Fig. [5l It is seen that for ordinary 
QNMs, /c(’’*^®)(s) and k'{s) practically coincide, whereas there are no resonances associated with 
hidden QNMs (#1,5,6). 


QNMs. The ordinary QNMs whose real parts of the complex-valued eigenfrequency, 
Rei^(mod) JR given frequency interval, can be determined from the transmittance 

spectrum T{u) of ID disordered samples, because each peak in the spectrum corresponds to 
a frequency whose value practically coincides with Reci;*^™°‘^h Moreover, when disorder 
is strong enough, so that L > ii^c, the distribution of the transmitted wave intensity along 
the sample reconstructs very closely the shape of the intensity of ordinary QNM eigenfunc¬ 
tions. In contrast, a hidden QNM is invisible (this explains the origin of the term ’’hidden”) 
in the transmittance spectrum and its intensity distribution is indistinguishable from that 
at a non-resonant frequency. 


Note that although the hidden modes are not displayed in the amplitude of the trans¬ 
mission coefficient, they are manifested in the phase of the transmission coefficient. The 
density of states at a frequency u is proportional to the derivative with respect to frequency 


of the phase of the complex transmission coefficient 


Our numerical calculations show 
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that each hidden mode adds tt to the total phase shift of the transmission coefficient exactly 
in the same way as ordinary QNMs. 



0.0 0.1 0.2 0.3 



Distance, a.u. 


FIG. 7. Difference between ordinary and hidden QNMs. (a) transmittance spectrum, (b) 
distribution of the incident wave intensity into the sample as a function of frequency and distance. 
There are two QNMs with nearby real parts of eigenfrequencies Rea;^“°^\ marked by dashed lines 
in the panels (a) and (b). Distributions of the intensities of eigenfunctions of hidden and ordinary 
QNMs along the sample are shown by thick blue lines in panels (c) and (d), correspondingly. 
Hatched red areas in panels (c) and (d) show intensity distributions of the incident waves whose 
frequencies coincide with of hidden and ordinary QNMs, corresondingly. 


The evolution of a hidden QNM as the degree of disorder grows is analogous to the 
evolution of a mode in a regular resonator when one of its edges becomes less transparent. 
This means that a hidden mode may be transformed into an ordinary (i.e., made visible in 
the transmission) by increasing the reflectivity of the corresponding edge of the sample, as 
illustrated in Fig. |8l 

The sample, whose transmission spectrum is shown in Fig. [5l contains three hidden 
QNMs 1,5, and 6) in the given spectral range. Distributions of the intensity I{n) {n is 
the layer number) for QNMs 7)^1 and 6 are similar to the distributions in resonators with 
right reflection coefficient smaller than the left reflection coefficient tl, < 1. The 

intensity distribution of QNM ^5 is characterized by the opposite inequality > 1. 

When the value of the right-end reflection coefficient is increased, new resonances appear 
in transmission for the initially hidden modes and 6, while mode #5 remains hidden 
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(Fig. [8^). In contrast, increasing the left-end reflection coefficient transforms QNM into 
a ordinary mode, whereas QNMs and 6 remain hidden in the transmission spectrum 
(Fig. [H)3). 



FIG. 8. Hidden modes 7^1,5, and 6, which are shown in Fig. [5l can be transformed in the ordinary 
ones by additional end reflections, (a) - additional right-end reflection transforms the hidden modes 
7^1 and 6 into the ordinary modes; mode remains invisible in the transmission spectrum, (b) 
- in contrast, the left-end reflection transforms mode ^^5 into a ordinary mode, while modes 
and 6 remain invisible. 

Important to stress that the separation of QNMs into two types, ordinary and hidden, 
occurs already at a very small disorder strength, s —)■ 0, when the localization length is 
larger than the sample length, £ioc L. 

The ensemble-averaged of the ratio of the number of transmission resonances, iV^es, which 
is the number of ordinary modes, to the total number of QNMs, iV^od, has been calculated 
numerically for a variety of randomly layered samples with different types of disorder (ran¬ 
dom reffection coefficients of the layers, r^, and/or random thicknesses dj, with rectangular 
and Gaussian distribution functions) in broad ranges of the disorder strength s and of the 
length of the realizations N. 

Figure |9] shows the average of Nj-^s/N^od over 10^ random realizations as a function of 
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the ensemble-averaged transmission coefficient (T) [panel (a)], and as a function of ratio 
of N to the localization length nioc (measured in numbers of layers), N/n\oc [panel (b)] for 
samples with N = 50,100,150, and 200 layers. At this scaling, all functions A'res((7"))/A^mod 
and A^res(A^/?^ioc)/A^mod for samples of different lengths merge in a single curve. 



(r) 



FIG. 9. Ratio Ares/A^mod as a function of the ensemble-averaged transmission coefficient (T) 
[panel(a)], and as a function of ratio of N to the localization length nioc [panel (b)] for systems of 
various lengths (number of layers N). The horizontal dashed red line marks y^2/5. 

It is seen in Fig. [9] that the difference between A^es and N-^oA appears when nioc S> A, 
and the ratio Ares/Amod varies weakly even when riioc A. Moreover, independently 
of the samples parameters, the average ratio Ares/Amod tends to the constant ■\/2/5 when 
TT-ioc oo. Thus, the existence of hidden modes and the universality of their relative number 
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is a general feature of ID disordered systems not specifically related to localization. 


IV. MEASUREMENTS OF TRANSMISSION EIGENCHANNELS AND TRANS¬ 
MISSION RESONANCES IN MULTICHANNEL SYSTEMS 


It is of interest to explore the ratio of the numbers of local maxima in transmission and 
QNMs in random multichannel systems and to compare to results for ID systems. We 
consider quasi-one dimensional (quasi-lD) samples with reflecting sides and transverse di¬ 
mensions W much smaller than the sample length L and localization length £ioc = -Whan^, 
hh < ^\oc■, L. Here, Whan is the number of channels or freely-propagating transverse modes 
in the perfectly conducting leads or empty waveguide leading to the sample and (. is the 
transport mean free path. The incident channels are thoroughly mixed by scattering within 
the sample. In contrast to transmission in ID samples with a single transmission chan¬ 
nel, transmission through quasi-lD samples is described by the field transmission matrix t 
with elements t^a between all Whan incident and outgoing channels, a and b, respectively. 
From the transmission matrix, we may distinguish three types of transmission variables in 
quasi-lD samples: the intensity Tba = the total transmission, and 

transmittance, T = J2at=i |tfea|^- The transmittance is an alog ous to the electronic conduc¬ 


tance in units of the quantum of conductance e^/h [u, 15, 40|. The ensemble average value 
of the transmittance T is equal to the dimensionless conductance, g = (T), which character¬ 
izes the crossover from diffusive to localized waves. In diffusive samples, the dimensionless 
conductance is equal to the Thouless number, g = 6 and the localization threshold is reached 
when g = 6 = 1 10, 12|. 

Signihcant differences between results in ID and quasi-lD geometries can be expected 
since propagation can be diffusive in quasi-lD samples with length greater than the mean 
free path but smaller than the localization length, i < L < Wc = Wvian ^, whereas a diffusive 


regime does not exist in ID since Wc = ^ 


4l| . For diffusive waves, QNMs overlap spectrally 


and may coalesce into a single peak in the transmittance spectrum. Thus we might expect 
that the QNMs within a typical linewidth form a single peak in transmission so that the 
ratio Wes/^mod is the ratio of the mode spacing to the mode linewidth. The mode linewidth 
is related to the correlation frequency in the transmission spectra, but the mode spacing 
cannot be readily ascertained once modes overlap. 
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The transmittance can also be expressed as T = '^^= 1 ’^ where the are the eigenval¬ 
ues of the matrix product tt^ 15| . The transmission matrix provides a basis for comparison 
between results for ID and quasi-lD, which is often more direct than a comparison based 
on QNMs, since the statistics of the contribution of different modes to transmission is not 
well-established, whereas the contribution of different channels is simply the sum of the 
transmission eigenvalues. In addition, transmission eigenchannels are orthogonal, whereas 
the waveform in transmission for spectrally-adjacent modes are strongly correlated 8j so 
that the transmission involves interference between modes. 

The transmission eigenvalue may be obtained from the singular value decomposition of the 
transmission matrix, t = UAV^ Here, U and V are unitary matrices and A is a diagonal 
matrix with elements yTy. The incident helds of the eigenchannels on the incident surface, 
Vn, which are the columns of V, in the singular-value decomposition are orthogonal, as are 
the corresponding outgoing eigenchannels, Un. Only a fraction of the Nrhan eigenchannels 
contribute appreciably to the transmissio n 114 1. In diffusive samples, the transmission is 
dominated by g channels with > 1/e jig . |43| . while a single eigenchannel dominates 
transmission for localized samples. The statistics of transmission depend directly on the 
participation number of transmission eigenhannels, M = '^n 2^1- Af is 


equal to 3g/2 2^ for diffusive waves and approaches unity in the localized limit 29|, [SOj 


Numerical simulations 


To explore the ratio A^res/A^mod over a broad range of g = (T) for multichannel disordered 
waveguides in the crossover from diffusive to localized waves, we carry out numerical simu¬ 
lations for a scalar wave propagating through a two dimensional disordered waveguide with 
reflecting sides and semi-inhnite leads. For diffusive samples in which there is considerable 
mode overlap since 6 = Su/Au > 1 12j, {Su and Az/ are the linewidth and the distance 
between spectral lines) the density of states (DOS), and from this the number of QNMs 


within the spectrum, can be obtained 
phase of the transmission eigenchannel 


rom the sum of the derivatives of the composite 


44l |. The derivative of the composite phase of the 


nth eigenchannel is equal to the dwell time of the photon within the sample in the eigen¬ 
channel. The total number of modes Amod in a given frequency interval is then the integral 
over this interval of the DOS. This has allowed us to determine the ratio Ares/Amod in the 
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crossover to localization. 

Simulations are carried out by discretizing the wave equation 

V'^E{x, y) + kle{x, y)E{x, y) = 0 


( 3 ) 


on a square grid and solved via the recursive Green function method Here, is the 
wave vector in the leads. Also, e{x,y) = 1 ± Se{x,y) is the spatially-varying dielectric 
function in the disordered region with 6e{x, y) chosen from a rectangular distribution and 
e = 1 in the empty leads. Reflections at the sample boundaries are minimal because the 
sample is index matched to its surroundings. The product of ko at 14.7 GHz and the grid 
spacing is set to unity. In the frequency range studied, the leads attached to the random 
waveguide support Achan = 16 channels which are the propagating waveguide modes. In 
our scalar quasi-lD simulations for a sample with a width W, the number of channels at 
frequencies above the cutoff frequency is the integer part of 2W/\. These channels should 
not be confused with the QNMs of the random medium which correspond to resonances of 
the medium with Lorentzian lines centered at distinct frequencies. In the simulations, the 
length of the sample L is equal to 500 in units of the grid spacing except for one deeply 
localized sample with g = 0.12, for which L = 800 and the width of the sample W is IGvr. 
Typical spectra of intensity, total transmission and transmittance are shown in Fig. [TOl for 
a diffusive sample with = 2.1 and for a localized sample with g = 0.3. 

We hnd that the numbers of peaks in the spectra of intensity, total transmission and 
transmittance in a single sample are nearly the same for each of the samples shown in 
Fig. [ini This is seen to be the case over a wide range of (T) in Fig. [TT] 

The DOS and so the number of QNMs within the spectrum in the samples of the same 
size are not affected by the strength of disorder so that the decreasing ratio N^es/Nmod 
with increasing (T) reflects only the decreasing number of peaks in the transmission spectra 
due to the broadening of the modes and the consequent increase in their spectral overlap. 
Since there are typically 6 QNMs within the mode linewidth for diffusive waves, we might 
expect the ratio N^es/Nmod to fall inversely with M, Ares/Amod 1/6 1/^ 3/2M. For 

deeply localized waves, however, this ratio is expected to approach unity as M approaches 
unity. This suggests that Areg/Amod ~ 1/M. in this limit. A plot of 1/M in Fig. [TT] 
shows that towards the diffusive and localized limits 1/M is close to the ratio Ares/Amod- 
For diffusive waves, the intensity correlation frequency does not change as the width of the 
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FIG. 10. Spectra of intensity, total transmission and transmittance for a localized sample drawn 
from a random ensemble with g = 0.3 (a)-(c) and a diffusive sample taken from an ensemble with 
g = 2.1 (d)-(f). Sharper spectral features are observed and spatial averaging is seen to be less 
effective in smoothing the spectra for localized waves than for diffusive waves. 


sample changes for fixed length and scattering strength since it is tied to the time of the 


flight distribution, which is independent of W 


46| . Since iVres is essentially the width of the 


spectrum divided by the correlation frequency of the intensity, the number of peaks within 
the intensity spectrum does not change. However, g and the DOS are proportional to iVchan, 
so that M increases with sample width and N^es/N^od is inversely proportional to M. In 


addition, the propaga 


approaching unity 


30| 


ion in a multichannel disordered sample is essentially ID, when M is 


These results suggest that a comparison can be made between propagation in both ID 
and multichannel systems via the ratio of the number of peaks in the transmission spectra to 
the number of modes normalized by M, N^^alM). This ratio may be expected to be 
close to unity for L S> ^loc- We consider the variation with g = (T) of the ratio MA^res/-^mod 
in quasi-lD and compare this with the corresponding ratio in ID in which M = 1. The 
values of this ratio in quasi-lD and ID are close, as seen in Fig. [121 
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FIG. 11. Variation of the Nj-es/N^o^ for transmission, total transmission and transmittance vs. 
g = (T) for multichannel random samples in simulations. The ratios obtained from microwave 
measurements of spectra of the three transmission variables in a multichannel localized sample 
with g = 0.37 are shown as the cross symbols and are in good agreement with the simulations with 
a similar value of g. The value of 1/M found in the simulations is shown as overturned triangles. 

B. Microwave experiment 

For quasi-lD samples in the crossover to localization in which spectral overlap is mod¬ 
erate, it is possible to analyze the measured held spectra to obtain the central frequencies 
of the QNMs and to compare these to peaks in transmission. Spectral measurements of the 
transmittance T were made in a copper waveguide of diameter 7.3 cm and of length 40 cm 
containing randomly positioned alumina spheres with index 3.14, over a random ensemble 
for which g = 0.37. The empty waveguide supports Nrh^r, ~ 30 propagation channels in the 
frequency range of the experiment: 10.0-10.24 GHz. The transmission matrix is determined 
from measurements of the held transmission coefficient between points on grids of 49 loca- 
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FIG. 12. Number of peaks in the transmission spectra per effective transmission eigenchannel, 
Nj-es/(N^ofi/Al), is plotted as a function of T = gjM. The quantity g/M is the effective transmis¬ 
sion coefficient per effective transmission eigenvalue of the quasi-lD system. Such normalization 
of the conductance g in quasi-lD samples makes possible a comparison with ID systems. The red 
line corresponds to a ID system; the experimental data is shown by the asterisk; the blue dots 
show the results of numerical simulations; and dashed line is drawn at the level ■\/2/5. Beyond the 
diffusive regime the ratio plotted rises towards unity for ballistic propagation. For ballistic waves, 
each of the ^channel transmission eigenvalue is unity so that the transmittance is N'channei and all 
eigenchannels contribute equally to the transmittance so that M = A'^channeb yielding g/M = 1. 


tions for the source antenna and detection antennas on the input and output surfaces of the 
waveguide for a single polarization with a grid spacing of 9 mm. Such measurements of the 
transmission matrix in real space for a single polarization are incomplete. The distribution 
of transmission eigenvalues determined from these measurements may differ from theoretical 


calculations 


42 


47j . We find, however, that the impact of incompleteness upon the statistics 


of transmittance and transmission eigenvalues is small as long as the number of measured 
channels is much grea ter than M, as is the case in these measurements of transmission in 
localized samples 30(]. Here M = 1.23 and therefore the statistics of transmission are not 
affected by the incompleteness of the measurement SOj. The influence of absorption in these 
samples is statistically removed by compensating for the enhanced decay of the held due 
to absorption Different random sample conhgurations are obtained by briehy rotating 
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and vibrating the sample tnbe. The probability distribution of the transmittance is in good 


agreement with the distribution calculated for this value of g 


30 
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m- 


We hnd the central frequencies and linewidths of the QNMs within the frequency range 
of the measurements by carrying out a modal decomposition of the transmitted held. A 
given polarization component of the held can be expressed as a sum of the contributions 
from each of the QNMs: 


E{r,uj) 




Tm/2 


^m/2 T i (cj ^m) 


( 4 ) 


Here am(r) are comp lex-valued amplitudes of QNMs. 

The central frequencies Um and linewidths T^ of the modes are found by simultaneously 
htting 45 held spectra. The transmittance as well as the Lorentzian lines for each QNM 
normalized to unity and the DOS, which is the sum of such Lorentzian lines over all QNMs 
are shown in Fig. [T3]for a single random conhguration. The DOS curves for diherent modes 
are plotted in diherent colors so that they can be distinguished more clearly. The DOS is 
also determined from the sum of the spectral derivatives of the composite phase of each 
transmission eigenchannel and plotted in Fig. [131 The DOS determined from analyses of 
the QNMs and of the transmission eigenchannels are seen to be in agreement. The dashed 
vertical lines in Fig. [13] are drawn from the peaks in the transmittance spectra in (a) to the 
frequency axis in (b). As found in ID simulations, each peak in T is close to the frequency 
of a QNM, but many QNMs do not correspond to a distinct peak in the transmittance. 
Frequently, more than one QNM falls within a single peak in T. 

The ratio of the number of peaks in spectra of transmittance to the number of QNMs 
averaged over a random ensemble of 40 conhgurations is 0.61, with a standard deviation of 
0.057. This is indicated by the cross in Fig. [TT] and is consistent with values of the ratio 
found in computer simulations. This value of this ratio is slightly smaller than the value 
0.65 found in simulations for ID sample with (T) = 0.37 found in ID simulations, as seen in 
Fig. [9l This may be attributed to the value M = 1.23 being larger than the value of unity 
in ID. This reflects the tendency of the ratio to decrease with increasing M as found for 
diffusive waves. 

Equation (|1|) offers an explanation for the fact that the number of transmission resonances 
can be smaller than that of all QNMs. If, for example, the transmitted held is a sum of 
two modes, from Eq. (|1|) it follows that the number of peaks in the transmission spectrum 
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FIG. 13. (a) Spectrum of transmittance T and the individual modes. The integration of each 
Lorentzian curve in the lower panel over the frequency yields the density of state of unity. There 
are 22 local maxima in the spectrum of T and the number of modes are 39. (b) Spectrum of 
the density of states. The sum of all the Lorentzian curves above gives the DOS of the sample, 
which is seen to be in good agreement with the DOS [panel (b)] obtained via the summation of 
the composite phase derivatives of each transmission eigenchannel. 

is either one or two, depending on the widths of the modes. 

C. Spatial intensity distribution of QNMs within quasi-lD disordered samples 

In order to fully characterize the QNMs and their relationship to peaks in transmittance 
in quasi-ID samples, it would be desirable to examine the longitudinal prohle of QNMs 
within the media. Because we do not have access to the interior of the multichannel sample, 
however, we explore the spatial prohle of QNMs using numerical simulations based on the 
recursive Greens function technique. The Greens function between points on the incident 
plane tq and within the sample r' can be expressed in a manner similar to Eq. (jl|) as a sum of 
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contribution from each of the modes, We hnd in the simulations that the spatial distribution 
of the mth mode obtained by decomposing the held into QNMs depends weakly upon the 
excitation point tq. We therefore average the spatial prohle for each QNM over the prohles 
obtained for all excitation points on the input of the sample. 

We consider propagation in a sample drawn from an ensemble with a value of g which is 
below unity but still not too small. In this case, QNMs still overlap but it is yet possible to 
analyze the held into QNMs. We present in Fig. [TUthat a spectrum of transmittance in a 
sample conhguration chosen from an ensemble with g = 0.26 and (M) = 1.16, together with 
prohles of a ordinary and a hidden mode within the spectrum. The nature of propagation in 
the sample might not diher appreciably from propagation in ID samples, for which M = 1. 
We hnd that the intensity distributions integrated over the transverse direction of the hidden 
mode in the transmission spectrum of the quasi-lD samples falls monotonically within the 
sample, while the ordinary mode associated with peaks in transmission is peaked in the 
middle of the sample. 

V. ANALYTICAL CALCULATIONS OF Wes/Y^od 

To calculate the average number of TRs in the limit s 1, we use the single-scattering 
approximation and write the total rehection coefficient r{k) of a ID system as: 

r{k) = J:^^^rnexp{2ikxn), (5) 

where Xn is the coordinate of the n-th scatterer. The values fcmax, at which the transmission 
coefficients, T{k) = 1 — |r(/c)p, has a local extrema, are dehned as the zeros of the function 
f{k) = d\r{k)\ydk = 2Re [r{k)dr*{k)/dk]: 

fikrn..) = = 0 . ( 6 ) 

Assuming hrst that ddi = 0, we obtain 

f{k) oc - n) sin [2k{m - n)do] 

= sin {2kldo) 

+T.^=irn-irj} = E/1i sin (2/c/do) ai- (7) 

Equation ([7j) is the trigonometric sum sin (uik) with “frequencies” z/; = 2/do and 

random coefficients ai. The statistics of the zeroes of random polynomials have been studied 
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FIG. 14. Transmittance spectrum and intensity distribution of QNMs in QID disordered samples, 
(a) Spectrum of transmittance T for a localized sample drawn from an ensemble with g = 0.26. 
The green circles indicate the central frequencies of the QNMs found from a modal decomposition 
of the field. Longitudinal intensity distributions of two modes with central frequencies indicated 
by the arrows in (a) are shown in (b) and (c) on a semi-log scale. The average intensity shown 
is integrated over the transverse dimension of the sample. The spatial profile of the mode with a 
peak in transmittance is seen to have a peak in the average intensity in the interior of the sample, 
while the intensity falls into the sample for a hidden mode. This resembles the behavior of hidden 
modes found in ID samples shown in Fig. [71 


in 


28l |. where it is shown that the statistically-averaged number of real roots Nj-oot of such 


sum at a certain interval Ak is 


Ak 

where af = Var(az) is the variance of the coefficients = 
the refiection coefficients are uncorrelated, then 


( 8 ) 


'L^^lrn+irnl + '^n=i^n-irnl- When 


Var(a;) = 2{N — 1)1^ (cTq -|- , 


(9) 


where cTq = Var(r) and r is the mean value of rj. The sums in Eq. ([8]) can be calculated 
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using Eq. ([9]), which yields [52|: 


= 2 (cT^ + 2rV^) - /) ~ - (a^ + 2rV^) 

(a^ + 2fV2) /^(A^ - /) ~ + 2r^^o) • 


( 10 ) 


From Eqs. (jH]) and (ITOi) we obtain 


2AkNdo 2 AkL 2 

iV„o. = ^^\/j = 2 —\/-. 


( 11 ) 


where L = Ndo- Since the number of minima of the reflection coefficient is equal to the 
number of TRs, A^res = -^root/2, and the number A^mod of QNMs in the same interval Ak is 
^mod = A/cL/tt, from Eq. flTT]) it follows that 


Ares/^mod — ^2/5. 


( 12 ) 


Although this relation was derived for systems with random reflection coefficients and 
constant distances between the scatterers, it also holds for samples in which these distances 
are random {6di 7 ^ 0). In this case, the frequencies u = 2ldd in Eq. ([^ should be replaced 
hy p = 2\xm — XTn±i\- Since the main contribution to the sums in Eq. ([H]) is given by the 
terms with large I ~ N, the mean value of \xm — Xm±i\ can be replaced by /do, in fhe case 
of a homogeneous distribution of the distances dn along the system. This ultimately leads 
to the same result Eq. flT^ . 


VI. HIDDEN MODES: SIMPLE MODEL 


In Sec. II, QNMs were introduced as solutions of the wave equation satisfying the outgoing 
boundary conditions. Their eigenvalues of QNMs, = k' — ik", can be calculated as 

roots of the equation M 22 = 0, where M is the transfer matrix, which connects waves’ 
amplitudes at the left and right sides of the whole system. The transfer matrix of the 
system which consists of A^ + 1 scatterers separated by N intervals has the form: 

M = • • • S^f^S.f,. (13) 


Here 


^ikdi 

0 

0 

Q-ikdi 


(14) 
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and Tj is the transfer matrix through the ith scatterer. Assuming that reflection and trans¬ 
mission coefficients are real, T* can be presented as 


1/ti -TijU 
-Tijti l/ti 


1 

u 




where I is the unit matrix, and (T 2 is the Pauli matrix. 

Omitting denominator matrix M can be written as ordered product 


M ^ P[ (/ - sri&^ Si, 


( 15 ) 


( 16 ) 


where substitution r* svi is used. Eq. (1T6|) allows presenting transfer matrix as a power 
series in s 1: 

N Af+l 

M = 1[S, + J2 5 "^, ( 17 ) 

71=1 71=1 

where matrix An contain various ordered products of matrices Si and n Pauli matrices. 
Matrices Si are diagonal, whereas Pauli matrix is anti-diagonal, so only even combinations 
of Pauli matrices contribute to M 22 - Thus, 

[(V+l)/2] 

M22 = e-^^^+ Y. (18) 

771=1 

where L is the sample length, and the coefficients contain various combinations of prod¬ 
ucts of 2m reflection coefficients r* with exponential multipliers exp(—ifcL -|- 2ikdij), where 
dij are the distances between any ordered pairs of scatterers. 

Neglecting terms with higher than powers in Eq. flT8|) . the dispersion equation, which 
dehnes eigenvalues k, can be presented as follows: 






CijC 


2k"dj 


= - 1 , 


( 19 ) 


where the coefficients Cij = are formed by various pairs of the scatterers, as it is 

schematically shown in Fig. [TSh . The greater is the distance dij between the scatterers, the 
larger are exponents in Eq. ([19]). 

The largest exponents are associated with the pairs of scatterers placed near the opposite 
ends of the sample. When k''do 1, there are many such pairs, located in blue regions 
in Fig. [T^ . whose associated exponents are of the same order of magnitude, . Let us 
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(a) 



FIG. 15. (a) - coefficients Cij are formed by various pairs of scatterers; (b) - coefficient ci contains 
all possible pairs of scatterers, linked schematically by red line ^1, coefficient C 2 is formed by the 
scatterers from blue and green regions, connected by red lines ^2. 

combine all such pairs in Eq. fll9p in one term ci and characterize them by one common 
exponent The number of scatterers near the sample ends, which form this group, can 

be estimated as Ueff — 3> 1, so that the lengths of blue regions in Fig. [T5b are 

~ Ueffdo < L. 

The next group, C 2 , which is associated with the exponent of the order of i^-'^eado)^ 

consists of pairs of scatterers, one from green and another from blue regions in Fig. ITHb . In 
such a way, Eq. flT^ can be approximately presented as 

+ c2e^k"{L-n,sdo) + _ _ ( 20 ) 

Strictly speaking, the phenomenologically introduced number Uefr varies from group to group, 
but when ^ is large enough and it is possible to neglect its variation. 

The coefficients c„ in Eq. (|20|) are the sums of random vectors in complex plane. For 
any given sample the lengths of these vectors are hxed, whereas the phases varies from mode 
to mode, so that the magnitudes of the coefficients c„, been averaged over many modes, can 
be estimates as 

(|cn|) ^ \/(r’2)%efr ^ (21) 

Using Eqs. fl20ll and (EIJ), one can calculate value of k", averaged over many modes. 
When <C 1, {k"{s)) is large and the second term in the parentheses in Eq. (|20D is small 
as compared with the hrst one 1) and can be omitted. Then, the average 
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solution {k"{s)) of Eq. (120|) is 


{k”{s)) 


= — In 




^ 1 1 1 
- In rieff ~ — In 


s 


( 22 ) 


agrees well with the result of numerical 


2L 2L 

The dependence {k"{s)) described by Eq. 
simulations, presented in Fig. |3]by red line. 

Expression Eq. (I22ll describes averaged over many modes dependence k"{s), but for any 
given mode this dependence can be different. Indeed, since Ues N (for example, Uefr — 10 
for s = 10“® in the numerical simulation presented in Figs. [2] and [3])) fluctuation of the 
values of |c„| for different eigenmodes can be rather large. In particular, |ci| for a certain 
mode can be much smaller, than |c 2 |. Presenting Eq. (|20D in the form 

s „-2fc"L 

1 12g-2fc"neffdo ^ 

Cl 

has solution 

1 


CiS^ 


(23) 


it is easy to see that Eq. 


k" ~ 


2ue//do 


In 


(24) 


when exceeds some critical value 5^,^, 


2 2 l ~ |-1 

s 3>SCTit = |ci| exp 


In 


(25) 


^effdo 

Solution Eq. (I24l) is independent of s and represents the hidden QNM (see Fig. [3]). 

Recall that c„ are formed by different groups of the reflection coefficients. In general, 
the similar, independent of s, solutions of the dispersion equation appear when magnitudes 
some hrst coefficients c„ in Eq. fl20l) are small as compared with magnitudes of the next 
coefficients. 

In order to demonstrate that independent on s solutions of the dispersion equation indeed 
correspond to the hidden modes, let us consider the system composed of three scatterers 
only. The dispersion equation Eq. fl20|) for this system is 




(26) 


When all r, are of the same order of magnitude, r, ~ r, and s is so small [k"{s) is so large] 
that exp[/i;"(s)(ii^ 2 ] ^ 1, the solution k” of Eq. fl26|) is 


k" ~ 


In 


2{di + d2) s‘^r‘^ 


(27) 














Equations fl26|) and fl27|) are particular cases of the general formulas (120|) and (12^ . 


If, for example, |ri| is small as compared with |r 2 , 3 |, but 2> s^rit 
there is another solution of Eq. fl26|) : 


Irirgl ^exp 


r di+d2 j 

£2 

[ di 

ri 


k” 


1 , r2 
—- In — 
2di ri 


(28) 


This solution is independent on s, similarly to the solution Eq. fl2T|) . 

Result of numerical solution of Eq. (j26l) is shown in Fig. [161 Fig. IT7I demonstrates relation 



FIG. 16. Dependence k"{s). Vertical and horizontal red lines mark Scrit and k" defined by Eq. (j28p 


between real part of the QNMs’ eigenvalues k'{s) and position of the peak in the 

transmission spectrum. Note that hidden modes are invisible in the transmission spectrum 
even when s Scrit- 
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FIG. 17. k'{s) (black dots) and kres{s) (red dots). When s <C Scrit, the value of k" is the same for 
hidden and ordinary modes. Nevertheless, hidden modes are invisible in the transmission spectrum. 


VII. SUPERRADIANCE AND RESONANCE TRAPPING IN ID RANDOM SYS¬ 
TEMS 


The model introduced in the previous Section can be used to study the segregation of 
superradiant states and trapped modes in regular quantum-mechanical and wave structures 
and to illuminate the analogy between this phenomenon and existence of two types of QNMs 
(hidden and ordinary) in disordered systems considered above. Behavior of modes in regular 
open structures as the coupling to an environment is altered, has been intensively studied in 
condensed matter physics, optics, and nuclear, atomic, and microwave physics. Common to 
all these studies is the appearance of two time scales when the coupling to the environmem 


via open decay channels increases and QNMs begin to overlap 


53 


; for a review, see 




and references therein. When the coupling to the environment is weak, the lifetimes of all 
states tend to decrease as the coupling increases. As the coupling reaches a critical value, a 
restructuring of the spectrum of QNMs occurs leading to segregation of the imaginary parts 
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of the complex eigenvalues and of the decay widths. The states separate into short-lived 
(superradiant) and long-lived (trapped) states. This phenomenon is general and, by analogy 


to quantum optics 


60| and atomic physics 6ll-l63| , is known as the superradiance transition. 


In more complicated structures, such of those consisting of two coupled oscillating subsys- 


es, the superradiance 


57l | that strongly couple 


terns, one with a low and the other with a much higher density of sta' 
transition is closely related to the existence of doorway states 
to short-lived QNMs with external decay channels. 

It is important to stress that along with the pronounced similarities between the resonance 
trapping in many-particle quantum systems, open microwave cavities, etc., and between 
the “bidding” of some of quasi-normal modes in disordered samples there are substantial 
differences as well. In particular, resonance trapping happens in regular systems considered 
in 
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59| when the coupling of the large number of QNMs to a much smaller number 


of common decay channels increases. Without disorder, the samples that we consider are 
perfectly coupled to the environment (total transmission at all frequencies). Finite coupling 
appears due to disorder, as the result of the interference of multiply-scattered random helds, 
and the role of the coupling parameter is played by the strength of the scattering inside the 
system. 

To reproduce the superradince phenomena in disordered structures we modify the model 
slightly by placing the random sample between two reflectors with reflection coefficients tl 
and tr, located at distances Sl and Sr from the edge scatterers. For simplicity, we assume 
that Sr = Sr = S. These reflectors can be included in the dispersion equation Eq. (]20|) as 
additional scatterers as follows: 


= -1. (29) 


Here cq oc vr^r contains the products tr^rTi with corresponding exponential multipliers, the 
largest of which, exp [2/c"(L -|- 5)], is separated in the explicit form. 

To approach the conditions at which superradiance and resonance trapping occur, we 
consider below (in contrast to the previous sections) the evolution of the eigenvalues of a 
given sample with fixed s when tr^r —)■ 0. 

When the product ItrTrI is large, the hrst term in the parentheses dominates and the 
solution of Eq. fl29|) is 


A:" 


-In 

2{L + 2S) 


1 

rRTR 


(30) 


31 





















If (5 = 0, the magnitudes of the exponents in the first three terms are equal. When —)■ 

0, the magnitudes of the additional two terms decrease and the solutions of Eq. fl29|) tend 
to their solutions in the original sample (without end reflectors), as shown in Fig. [TSl 




FIG. 18. Two reflectors are placed at the sample ends, = tr = rend- Modes marked by letters 
correspond to the same modes in Figs. [21 [3l (a) - Trajectories of eigenvalues as the coupling grows, 
(b) - k"{rend)- The life time of the hidden QNMs decreases much faster than the life time of the 
ordinary ones. 
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When (5 7 ^ 0, the trajectories of the eigenvalues in the complex plane are more complicate. 
Although most of the eigenvalues finally reach the same positions as in the original sample, 
there are eigenvalues, for which k" —)■ oo as r^ ji —)■ 0 (see Fig. Indeed, the hrst two 
terms in Eq. always dominate when 5 7 ^ 0 and k" ^ 00 . In this case Eq. can be 
written as 

+ s {vlTn+i + rRVi) ~ 0, (31) 


where the largest term in cq, which corresponds to the largest distance L + 6 between the 
end reflectors and the sample scatterers, is explicitly presented. Solution of Eq. f^T]) 


k" 



rN+i n 

tr Tl 


(32) 


tends to infinity, when even one of the reflection coefficients tl^r —)■ 0 . 

The reason why these solutions “run away” when the coupling to the environment is 
maximal, is very simple. The original system without end reflectors has AfcL/vr eigenmodes 
in the given interval Ak, whereas the same system surrounded by the reflectors has Ak{L + 
25)/t: eigenmodes in the same interval. Thus, some of modes should leave this interval Ak 
when the system returns to its original state. 

The superradiant transition in periodic and disordered quantum system, which consist 


of a sets of potential wells, was studied in 


53| using effective Hamiltonian approach. It was 


shown that the transition occurs when the coefficient 7 , which characterizes the coupling 
with an environment, reaches the value of the coupling H between the wells, 7 ~ H. In 
the considered above system 7 ~ 1 — \rR^i,\^ and H ~ 1 — so that the superradiant 

transition occurs when {tr^rI ~ This condition agrees well with presented in Fig.fT^ 

results. 

Hidden modes can be associated with superradiant states, while normal modes are 
trapped resonances. Thus, \/2jh ~ 0.63 correspond to the fraction of the modes which 
are trapped. This result agrees with where this value was estimated as > 0.58, and 
1 — 7/275 ~ 0.37 is the fraction of the modes which are superradiant. Note, that the original 
disordered sample is already coupled to the environment, so that the coupling strength is 
limited by the intrinsic properties of the sample and cannot exceed this value, even when 
the end reflectors are fully transparent. 
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FIG. 19. The same as Fig. [T^ but reflectors are placed at a some distance from the sample ends. 
There is one eigenvalue whose k" grows unlimitedly 


VIII. CONCLUSIONS 

In conclusion, we have studied the relationship between spectra of quasi-normal modes 
and transmission resonances in open ID and quasi-lD systems. We start from homogeneous 
samples, in which each TR is associated with a QNM, and vice versa. As soon as an 
arbitrarily weak disorder is introduced, this correspondence breaks down: a fraction of the 
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eigenstates becomes hidden, in the sense that the corresponding resonances in transmission 
disappear. The evolution of the imaginary parts of the eigenfrequencies of the hidden QNMs 
with changing disorder is also rather unusual. Whereas increasing disorder leads to stronger 
localization of ordinary modes so that their eigenfrequencies approach the real axis, the 
imaginary parts of the eigenfrequency of hidden modes changes very slowly (and may even 
increase when external reflectors are added to the edges) with increasing disorder, and 
begin to go down only when the disorder becomes strong enough. For weak disorder, the 
averaged ratio of the number of transmission peaks to the total number of QNMs in a 
given frequency interval is independent of the type of disorder and deviates only slightly 
from a constant, ■\/2/5, as the strength of disorder and/or the length of the random sample 
increase over a wide range. This constant coincides with the value of the ratio Wes/Nmod 
analytically calculated in the weak single-scattering approximation. As the strength s of 
disorder keeps growing, ultimately all hidden quasimodes become ordinary. This means that 
in ID random systems there exists a pre-localized regime, in which only a fraction of the 
QNMs are long-lived and provide resonant transmission. If the coupling to the environment 
is tuned by external edge reflectors, the superradiace transition can be reproduced. In 
quasi-lD samples, a genuine diffusive regime exists in which some QNMs coalesce to form 
a single peak in transmission with width comparable to the typical modal linewidth. In 
such samples, hidden modes have been discovered experimentally and their proportion of all 
QNMs in the crossover from diffusion to localization was fairly close to the same constant. 
The number of peaks in spectra of transmission, as well as in total transmission and in 
transmittance are nearly the same and fall well below the number of QNMs. Though the 
ratio Wes/A^mod may be small, we hnd in microwave experiments and numerical simulations 
that once the number of QNMs is divided by the effective number of channels contributing 
to transmission to give MNp/Nm, this function is similar to results in ID samples. 
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